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In various unified extensions of the Minimal Supersymmetric Standard 



Model, the Yukawa couplings of the third generation are predicted to be of 



the same order. As a result, low energy measured mass ratios, require large 
ratios of the standard model Higgs vacuum expectation values, correspond- 
ing to a large value of the parameter tanf3. We present analytic solutions for 
5_i ■ the Yukawa couplings and the Higgs and third generation squark masses, 

in the case of large top and bottom Yukawa couplings. We examine re- 
gions of these Yukawas which give predictions for the top mass compatible 
with the present experimentally determined top mass and provide useful 
approximate formulae for the scalars. We discuss the implications on the 
Radiative Symmetry Breaking mechanism and derive constraints on the 
undetermined initial conditions of the scalars. 
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1 Introduction 



If one adopts the idea of unification of all forces, then the most natural candidate 
which is at the same time the simplest extension of the standard model (SM) 
of strong and electroweak interactions, is the minimal supersymmetric standard 
model (MSSM). Indeed, the low energy measured values of the three gauge cou- 
plings are consistent with simple unification values at an energy of O(10 16 GeV), 
provided that the SM - spectrum is extended to that of the MSSM at an en- 
ergy scale M e s {j SY ~ ITeV. Furthermore, unified supersymmetric models solve 
successfully the hierarchy problem [Ij of their corresponding non-supersymmetric 
versions. Recently, in view of the possible observation of supersymmetric particles 
in the upcoming colliders, there has been a growing interest in the determination 
of the complete mass spectrum of MSSM. Nevertheless, although the MSSM is 
considered as a successful extension of the SM, it still has a number of arbitrary 
parameters which are not fixed by the theory. It is expected however that these 
parameters could be reduced by embedding the MSSM in a more fundamental 
theory such as supergravity or superstrings. In the present status of MSSM and 
its Grand Unified extensions, (SUSY-GUTs), the Yukawa couplings and the ini- 
tial conditions for the scalar mass terms are treated as arbitrary parameters of the 
theory. Most of the SUSY-GUTs predict certain relations between the Yukawa 
couplings of the third generation, up to an overall constant which is left arbitrary. 
Certain string constructions may relate this constant to the unified value of the 
common gauge coupling at the string scale, usually up to a factor of 0(1). Most 
common relations that hold in these models are 

ht ~ hf, ~ h T , h t > hb ~ h T or h t ~ hi, > h T (1) 

Remarkably, it was shown by numerical analyses 0, that the top mass is predicted 
to be rat ~ 180GeV which is compatible with the recent experimental findings^]. 
Furthermore, in the context of MSSM, the top mass has been calculated in the 
recent literature mainly by numerical methods or semi-analytically in the case of 
the exact equality h t = hb |§, |6|. In the above cases, the experimentally observed 
low energy hierarchy m t 3> mj, demands a large ratio of the two Higgs vacuum 
expectation values (vevs), tan(3 = ^ 3> 1. 

As far as the scalar masses are concerned, the aforementioned theories have 
not been very predictive yet. In the simplest treatments, one usually assumes that 
all scalar masses (i.e. higgses and soft susy masses ) have the same boundary con- 
dition at the unification scale (universality), while their splittings at low energies 
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arise mainly from the renormalisation group running. In general, however, one 
expects that supergravity or superstrings will induce non-universal conditions for 
the scalar masses. 

In a certain class of supergravity models 0, each scalar mass rhi can be 
parametrised in terms of a coefficient & ~ C(l) (calculable in specific models), 
and the supersymmetry breaking scale which is proportional to the gravitino 
mass rrt.3/2 corresponding to an approximately flat direction of the underlying 
fundamental theory. In the effective low energy theory, the over all scale - set 
by the gravitino mass m 3 / 2 - is not known, since m 3 / 2 depends on singlet scalar 
field vevs of the higher unknown theory. It has been claimed however, that since 
m 3/ / 2 is a field dependent quantity, its value can be determined dynamically by 
minimizing || the vacuum energy of the effective low energy theory with respect 
to 777.3/2. This procedure has shown that under general assumptions for the higher 
theory and a wide choice of £i's||, the m 3 / 2 scale is of the required order, i.e. 
~ 0(mz). This tells us that supersymmetry is broken at the right scale to protect 
scalars from large radiative corrections, while the resulting susy scalar spectrum 
is accessible by the future colliders. 

While non-universality is not an inescapable prediction of the fundamental 
theory, it is a necessity for large tan/3. Indeed, if we are to take seriously pre- 
dictions of exact or approximate h t , hb equality at Mu, the radiative breaking of 
the SU(2) x U(l) gauge symmetry |L0| cannot be realized naturally if we assume 
universal boundary conditions - in particular for the higgs mass parameters. 

In the present work, we wish to explore the region of the parameter space 
of £j = m^/r?73/ 2 ,... which naturally leads to radiative symmetry breaking of 
the electroweak gauge group, protecting simultaneously the scalar quarks from 
receiving negative masses. This latter phenomenon is quite possible due to large 
negative radiative corrections on the scalars of the third generation. We should 
also point out that another possible consequence of using non-universal boundary 
conditions for scalars [D] is the appearance of enhanced flavor changing neutral 
currents (FCNCs) ||12|| . To suppress FCNCs within the experimentally accepted 
limits, it is adequate to consider only cases diagonal in the flavor space. In our 
analysis we are going to derive semi-analytic formulae which will prove useful as 
will be explained in the next section. 

Our paper is organized as follows: In section 2, in order for the paper to be self 
contained, we start with a short review of the low energy effective supersymmetric 
theory and explain our notation. We proceed by presenting the analytic solution 
for large h t — ht couplings and comment on their applicability to the minimization 
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conditions of the low energy potential with respect to higgs vevs and the gravitino 
mass. In section 3 we apply our solutions to derive analytic expressions for various 
useful quantities (including tan/3). Section 4 deals with the coupled differential 
equations for scalar mass of the third generation in the regime of large tan/3. It is 
shown that the solution of the system of equations for all third generation scalar 
masses is reduced to one simple first order equation of Riccati type. We derive 
analytic solutions in the limiting cases of low and large tanf3 cases and explore 
the non-universal conditions which are compatible with RSB. Finally, in section 
5 we present our numerical results and conclusions. 
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2 Large Top-Bottom Couplings in MSSM. 



The theory which concerns us here is a supersymmetric low energy effective theory 
based on the general structure of spontaneously broken N = 1 supergravity |L3] . 



Although such a theory still leaves a lot of arbitrariness, in the minimal case the 
effective theory can be described in terms of five new parameters only, namely 
the gravitino mass term 777.3/2 the m^, \x , m 2 = Bfi mass parameters and the 
trilinear parameter A. In the present paper however, for reasons explained already 
in the introduction, we will assume non-universality for the higgs and the other 
scalar mass parameters at the unification scale Mjj, thus our parameter space 
will be enlarged. At the tree level, the masses of the fermions are obtained from 
the trilinear superpotential given by 

W = I o'j U -i + h%Qid)H x + h i lL l e)E x + ^H X H 2 (2) 

Supersymmetry is broken softly through the " soft - breaking " terms which 
provide with masses all the s-particles. Thus, the gauginos A» receive masses 
from terms -mjAiAj. Scalar mass terms arise from the effective potential which 
in the minimal case is of the form 

£ = S|^| 2 + Am 3 / 2 (W + W*)+m2 /2 S| 0, | 2 + Bm^^—- + h.c.) + ... (3) 

where the {...} represent D-term contributions. A, B are parameters of 0(1) and 
in the case of non-universality they may have different values for each superpo- 
tential term, i.e. A = {Aq, Ap, ...} etc. Thus, generalizing the Lagrangian in 
(0) for the non-universal case, one obtains soft squark and slepton mass terms of 
the form fh\q* n q n ^ where q n = {Q n , u c n1 d c n1 L n , e c n1 is n }. Moreover, due to the trilin- 
ear superpotential couplings in (||) there correspond soft-susy breaking trilinear 
couplings of the form 

A Q h%Q iU °H 2 + Anh^QidjHt + Al^L^H, + B/j l H 1 H 2 (4) 

where now A n = (A n + 3B n )m 3 / 2 , with n = Q,D,... and B = (A^ + 2B^)m 3 / 2 . 

The RG-improved effective potential for the neutral particles can be written 
as follows 



Vq(Q) = {m 2 Hl + ^ 2 ) | H x | 2 +{m 2 Ha + /i 2 ) | H 2 | 2 +m 2 (# 1 # 2 + h.c. 

h/2 



2 / 

+ 9 —-^-{\ H x | 2 - | H 2 | 2 ) 2 + m - V (Q)mt /2 (5) 
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where 



is the one loop contribution, necessary to stabilize the minimization with respect 
to the higgs vevs =< Hi >, while 

StrM 2 (z,z) = Y,(-V 2Sn (2sn + l)ml(z,z) (7) 

n 

is the sum over all particles with spin s n and mass m n . 

The term n(Q)my 2 comes from the underlying supergravity theory and can 
be interpreted |5) as a contribution to the cosmological constant. Notice that 
terms of 0(m\, 2 Mp l ) are not included since such contributions would destroy the 
hierarchy m 3 / 2 <C Mp;. Specific examples with these desirable features may be 
found in four dimensional string models |T5[ . 

The quantum corrections to the above potential depend on the particle spec- 
trum masses which are computed with the use of R.G.E.s in terms of the initial 
values at the unification scale. A precise knowledge of the R.G. effects in the 
low energy theory is important not only from the phenomenological point of view 
(measurements of sin 2 8w, a em , top mass, e.t.c.) but for theoretical reasons as 
well. In particular the Higgs mass parameters and squark masses of the third 
generation play a very important role in the radiative breaking of the electroweak 
symmetry on the one hand and the preservation of the £77(3) confining group 
and the electric charge on the other. 

From phenomenological analyses, two limiting cases have been of particular 
interest. The first, is when the two higgs vev's are of the same order (tan/3 ~ 
(9(1)) while the difference between the top mass and the other fermion masses is 
attributed to a large top Yukawa coupling relative to the bottom, h t ^> h^. 

It is interesting that this case can easily be treated analytically both in the 
fermion and scalar mass sector at the one loop level |Tj| [I7|]. Since h t is the only 
non-negligible Yukawa, all the others can be ignored and the relevant R.G.E.s 
decouple, leading to first order differential equations which can be solved by 
standard methods. The same procedure can be used for the scalars including 
those of the third generation. 

The case with a large bottom Yukawa coupling (tan/3 ^> 1) is much more 
complicated. In a previous work [fTg|] , however, it was shown that if one neglects 
the small differences in the U(l)y factors, an analytic solution of the h t — 
coupled differential system can be obtained easily. In particular, in the absence 
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of a large tau Yukawa coupling, one can write the h t — hb equations at the 1-loop 
as follows 



h\ = ±_{6 h *+hl-G Q }ht 



with 



dt 1 

jhl = ^{ht + 6hl-G B }hl (9) 
Gq = Ec^ 2 , G B = j2J B g* (10) 

i=l i=l 

where d Q = {f,3,g} and 4 = {f,3,^} for S77(3), SU(2) and £7(1) respec- 
tively. 

In order to solve Eqs(|8|-|9]), we define the new parameters x,y through the 
relations hf = •JqX, h\ = JqU, with 

7j = e W {^- 2 J t t0 G Q (t')dt'} 

J 

= n(i-^(*-w) # (id 

and make the following transformation 

U = J^yff^^ dI= £-^ dt (12) 

where the parameter k = 4:X y /(x — Uq) 7 ^ depends on the initial conditions 
Xq = hf and y$ = /^ . Then, we can form a differential equation for the new 
variable u which can be cast in the form [TT3] 

^L = -\kTn (13) 

V 1 4" u 6 

Then, by simple integration it can be shown that the solution can be given in 
terms of hypergeometric functions, 2-^1 (a, b, c; z), i.e. 

uVW > F 4> "TTP 4 ) - ^ /10 2 iM~ ^, 4 — , ) = ^^o 6/5 /(t) (14) 

2 10 10 — u 2 10 10 —Mo 12 

with /(t) = ^2 ft lQ(t')dt' . Thus, for a given set £o,2/o of initial conditions at 
Mjj we can calculate the u,uj(t) values from equation (14) and determine the 
h t , hb from 

hi = l 2 Q x (15) 



hi = IqV (16) 

In ref similar solutions where proposed. A further study of the solutions by 
these authors showed that all h ti b > 1 initial values accumulate on a curve on 
the h t , hi, plane. In the region h t m h^, this curve is slightly deformed when the 
difference of the hypercharge in the two RG-equations is taken into account. A 
numerical fit has shown the relation || 

(ht - 0.015) 12 + (hi + 0.045) 12 = (-^=) 12 (17) 

The main deformation occurs around the diagonal, i.e. when h t ~ hj,. The top 
mass solutions presented in the subsequent section are slightly affected by this 
change. In any case, this uncertainty on the top mass is in fact overwhelmed by 
much bigger uncertainties arising from sparticle loop corrections on the bottom 
mass. A more detailed discussion of the latter will be presented in the next 
section. For most applications, (in particular when u 1), it is appropriate to 
expand the hypergeometric function and obtain a simplified form of the above 
solutions. They can be expressed in terms of a simplified form of the function 

L0(t), 

u{t) ~ uWTymmW ( 8) 

where the value of the hypergeometric function at u = Uq is 2-^1° ~ 1 for the 
limit of interest, i.e. for ht,o ~ Vo- Formula ( |18D is going to be used later in the 
calculations for the scalar masses. Note finally (for later use) the useful relation 
between x, y variables 

(^M 7 ={— ) 6 - (19) 

In particular, in certain supergravity theories the Yukawas are subject to vari- 
ous constraints which may be combined with the above formulae and result to 
interesting predictions || [19 . 



The analytic h t — solution neglects small differences arising from the differ- 
ent hypercharge assignment of the top and bottom quarks. This approximation 
is quite reasonable since other uncertainties (neglect of Yukawa couplings, one 
loop approximation e.t.c) of the same order are also introduced in the analytic 
solution. In any case, it is clear that the vastly different values of m t and mj 
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masses cannot be attributed to effects due to these small corrections unless a 
considerable fine tuning of the parameters occurs. 

The usefulness of the analytic expressions for h t — hb couplings in the large 
tand regime, is largely exhibited in the case of the scalar mass calculations. As 
has already been mentioned, the knowledge of the scalar mass parameters affected 
by the large Yukawa couplings at low energies, is very crucial and decisive for the 
particular supergravity or string model under investigation [pOfl . To be more 
concrete, let us assume that only h t , Yukawas are large. In this case the up 
and down squark masses of the third generation as well as the two higgs mass 
parameters of the standard model {m 2 Hi1 i = 1,2) doublets receive large negative 
contributions proportional to integrals involving h 2 ,hl- If <C h t , and h t is 
large enough, then m 2 H2 mass parameter turns to a negative value at low energies 
triggering the SU(2)l(&U(1)y breaking radiatively. On the other hand, m 2 Hi mass 
does not receive large negative corrections since coupling is small. This is of 
course desirable since (as it turns out in this case), the neutral Higgs potential 
is stable and the vev's of the Higgses are finite as expected. Besides it can be 
proved |17j that for the experimentally accepted top quark mass and in the largest 
part of the (mi/2,%/2) parameter space the squared masses of the top-squarks 
receiving negative contributions are not driven negative, protecting in this way 
the color and charge quantum numbers. 

The situation is much more complicated in the case of large bottom coupling 
too. Then the other Higgs mass m 2 Hl may also be driven negative at low ener- 
gies. In order to make the analysis more transparent, let us write the conditions 
obtained from minimization of the Higgs potential, with respect to v^s 

ml = -^(m 2 Hi +m 2 H2 +2^ 2 )sm2f3 (20) 
9 m«- — m% n tan 2 8 „ . , 

mi = 2 ta^ 1 - 2 "' (21) 

In the large tand scenario that we are interested here [tand ^> 1) we may approx- 
imate the above formulae so that m 2 z ~ —2 (m 2 H2 + fi 2 ) and — m\ tand ~ m 2 Hl + 
m 2 H + 2ji 2 > 0. Combined together, these equations lead to m 2 H — m 2 H > m 2 z 

2 

and fil = fi 2 + m 2 Hi > ^f- > with the conclusion that the stability of the po- 
tential is ensured only for a positive value of /if mass. Another compelling reason 
is also the protection of the color and charge quantum numbers. The above de- 
scription shows clearly the necessity of obtaining analytic solutions of the scalar 
mass parameters in the case that tand ^> 1. 

A final comment we would like to make concerning the usefulness of the pro- 
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posed semi - analytic procedure is also based on theoretical motivations. Indeed, 
as we have already mentioned, a particularly interesting idea is that the vacuum 
energy density could be minimized with respect to the gravitino mass m.3/2 in 
order to fix the mass of the latter - and therefore the supersymmetry breaking 
scale - dynamically The minimization of the vacuum energy with respect to 
the 777.3/2 parameter leads to a certain condition for the renormalised scalar mass 
parameters and Higgs mass terms at low energies, which reads || 

2^i(Q)+77(QK /2 = (22) 

while the coefficient rj(Q) obeys a renormalisation group equation of the form 

dr](Q) 1 



dlog(Q) 1287T 



: StrM 4 (23) 



The dynamical determination of the 7773/2 scale requires the knowledge of the 
evolved soft (mass) 2 -parameters involved in the above constraint. Thus, analytic 
knowledge of their dependence on the boundary conditions are very useful. In 
particular, in the case of non - universal boundary conditions of the scalars at Mjj, 
the parameter space at M v can be represented by a vector £ = (fan £ff 2 > £qj •••) 
where mj(Mu) = £»m§/ 2 . It is not easy to handle numerically the above system 
in this case and therefore an approximate analytic solution of the scalar mass 
parameters would be very important. 

In what follows we wish to extend our previous results of the t — b analytic 
expressions. We give expressions for various useful quantities including tanf3 and 
derive simplified formulae for h t ^ Yukawa couplings. In particular, in the next 
section we develop the formalism which is going to be used in section 4 for the 
derivation of the scalar mass analytic expressions. 
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3 Phenomenological analysis of b — t solutions 



Let us define two new parameters a and r 

3 



eXP ^l/ idt '} (24) 
3 r* 



a = exp{— J hldH) (25) 
The ht-hb RGE's can be combined to give 

± = ^(lf (26) 
a y \(jJ 

where dots stand for derivatives with respect to t — ln(Q). Integration of the 
above gives 

111111 

= (27) 

Rewriting now h t — h b RGE's in terms of (cr,r) and substituting into (|2T| ) we 
obtain 

x-y 

TO - 



- y 

(28) 



6 

U \ 5 



u 

where uo = u(xo, yo). On the other hand, we can express the ratio ^ in terms of 
the integral 



~o = exp{^j\h 2 t -hl)dt} 



1 1 

(.r/>\ I (x-y)dl\ (29) 

*0 



which, with the help of C |T~3l , [T2j) , can be expressed as follows 



o \xoVTTTi-iJ 

The explicit dependence of a, r parameters on the new variable u defined in ( |i"2"l) 
is not only necessary for the scalar mass solutions we are going to discuss in the 
next section; it will also prove useful in other low energy parameters. Let us for 
example consider the evolution of the ratio of the two higgs vev's ( tan/3), in the 
case we are examining here, i.e. when only h t ^ couplings are large. In this case 
we have the following evolution equation 

d 13 

_ton/5 = -^^(h 2 t - h 2 b )tan(3 (31) 
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Integration of the above gives 



a 



tan{3 = tan(3 — 



, T 



Thus tan/3 evolution is also determined solely from the parameter u and the initial 
values of the two Yukawas. From (|32| ) it is obvious that for large u values the 
w-dependent term in ([32] ) tends to unity and the ratio Rp = tanf3 / tan(3 is given 
approximately by (^t,o/^b,o) 3 ^ 5 - 

In the most general case, it would be interesting to combine the above result 
with the (m t ,mf,) mass relation to eliminate the parameter u. Notice that to ob- 
tain this relation in the case of tan(3 ^> 1 one should take into account corrections 
to the bottom mass, arising from loop-graphs [f2~T| , |22 , p3fl containing supersym- 



metric particles. It is possible however to assume the existence of PQ- or R- type 
symmetries which may prevent the appearance of such large corrections [22]. For 



the sake of simplicity we will assume that due to some kind of symmetry such 
corrections are absent. In this case one finds 



h t ,o, n ftan(3 \ 
— tanon 



m t = m bT ^tan(3 ^ ^—^ j (33) 



h 



6,0 



where mb is the running bottom mass at the scale Q = m t , related to the bottom 
mass at Q ~ nib with a renormalisation group factor 7]b w 1.41, i.e., m^m;,) = 
r]i } mf ) (m t ) . In figure 1 we show the variation of rp = (tan(3 / tanPo) 2 ^ 3 vs the 
parameter u whose range (uo,uj) is determined from the initial values of the 
couplings h t fi,hbfi and the solution (JTJ). (uf corresponds to the value that u 
receives at the scale m t ). Since we are interested for large tan[3 values, we plot rp 
for several ^,0/^6,0 ~ ratios close to unity. In all cases we observe a small change 

2 

of Tp which approaches asymptotically its maximal value, (/fyo/^&.o)^ as u ~~ * 00 • 
Therefore, in these cases putting all together, one finds the simple relation 

m t ~ mb(mt)tan(3 (34) 

where mb(m t ) is the mass at the scale Q = m t . 

In Table I, we present selected values of the top mass, in the case where the top 
Yukawa coupling is close to its non-perturbative value, (h tt0 = 3.0) and various 
values of the t/b Yukawa ratio close to unity, h t /h b fi > 1, i-e. when ft, 6)0 Yukawa 
coupling is also large. Both, running and physical top masses are shown while in 
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the last column we present the corresponding values of tan ft, assuming a central 
value for the bottom quark, i.e. m b {m h ) = rib^hb cos ft = i.25GeV. Here, as 
already explained, n& includes the running from the scale Q ~ mt down to the 
scale rrib and is taken to be rib ~ 1.41. We start the running at the unification 
scale Mu, using as inputs the unification scale itself and the value of the common 
gauge coupling ay at M v together with the initial values of the Yukawa couplings 
ht,o, hb,o as shown in the Table I. Then, in obtaining the low energy values of a em , 
a 3 , and sin 2 9 w , we use the following ranges 

a em ~ l = 127.9 ± .1, a 3 = .12 ± .01, sin 2 9 w = .2319 ± 0.0004 

For the large h t values chosen here, the top mass is in remarkable agreement 
with the recent experimental data 0]. 

Of course the accuracy of the results of the table should not be overesti- 
mated. A more precise calculation should involve uncertainties arising from the 
strong coupling measurements and other parameters as the bottom corrections 
mentioned earlier. Our aim here is to show that our solutions result to stable 
m t , m t h , tan/3 predictions against h tt0 /hb t o variations. 



h t ,o 


ht,o/hbfl 


m t /GeV 


mf/GeV 


tanft 


3.0 


1.17 


176.7 


184.7 


58.22 




1.25 


177.0 


184.9 


58.10 




1.30 


177.1 


185.0 


58.06 




1.36 


177.2 


185.2 


57.90 




1.42 


177.3 


185.3 


57.78 



Table I. m t and tan/3 values for ht,o = 3.0 and various htfl/hbfl-ratios shown 
in the second column. We fix rrib = nb-^hbcosft = 4.25GeV . Large sparticle 
loop corrections on the bottom mass may change the above predictions. 

Indeed, as has been already noticed, uncertainties from sparticle loops contribut- 
ing to nib are not included. Therefore, the analytic solutions of the scalar masses 
(presented in the subsequent section) which are involved in the relevant loops 
may considerably simplify the examination of the above corrections. The Srrib 
corrections will have a direct impact on the determination of the angle ft and 



therefore on the top-mass. The relevant loop-graphs [21] lead to the following 



result m, H 

,2 
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where the function I(x, y, z) is given by 

xy\n(x/y) +yz\n(y/z) + xz\n(z/x) 

n x , V> z ) = 7 77 77 x , ( 36 ) 

[x - y){y - z)(x - z) 

In the above, m g is the gaugino mass, while m- bil mi i are the eigenmasses of 
sbottom and stop. A crucial role is played by the relative sign of the two terms 
involved in the above formula. The gaugino mass is given by its one loop for- 
mula, irtg = Our analytic procedure in the next section can also provide 
a simple formula for the trilinear coupling A t = f t (A tt o, Ab t o) — I' T m i/2> where f t is 
a linear function of A's while I' T an integral to be discussed in the subsequent sec- 
tion. For reasonable ^4(i,&),o initial values A t can be approximated A t « —I' T m\/2- 
Thus, there is a relation between the trilinear parameter and the gaugino mass, 
A t Pd —0.6rrig which can simplify our discussion. Therefore, there is a partial can- 
celation of the above two terms in the formula, but in the case tan/3 ~ (50 — 60) 
as the table shows, the correction is still large. 

A precise estimation depends on the details of the scalar mass spectrum in- 
volved in the integral functions I(x, y, z) and the specific initial conditions of the 
scalars. This is nevertheless possible for any set of initial conditions on the scalars 
with the help of the analytic solutions of the next section. Furthermore the sign 
of ^ plays also a decisive role. It can be easily seen that such corrections can 
lie in a wide range, from a few up to 40 %. Evidently, similar graphs to those 
considered above, can also result to corrections on the tau mass. Nevertheless, in 
this case it is found|23| that they are at most 5% and do not affect essentially the 
top mass prediction. It is worth noting here that in particular unified models as 
in the case of SU(5) x £/(l)P7|, the b — r Yukawa couplings are not necessarily 
equal at the Unification scale. Thus, if these corrections come with the right sign, 
they may reconcile large top mass values, with the condition hb > h T . In such 
cases, our present analytic procedure where the h T coupling is ignored in valid 
within all the parameter space. 
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4 Analytic Solutions for the 3 rd Generation Scalar 

Masses 

We have already pointed out that in unified supersymmetric models where only 
h t , h b couplings are large, the differential equations determining the evolution of 
the soft scalar masses of the third generation and the two Higgs mass parameters 
rriHi 2 are coupled. In the general case with non-zero trilinear parameters A t , At,, 
these equations take the form 



^-^Mli + ^MU^Ml + ^S (37) 

f^~E«»? + 4^-^ (38) 
1 Y,^Mh! + i$- 2 Ml + ^S (39) 



dm 2 Hr> 1 „ hi , j9 l«i 

dfh 2 D 1 ^ D„2 2, Ji , 1 «1 



dt - 8 ^c?MM + 2^MU~S (40) 



where 



^ = + mj} + m| 2 + A t 2 (42) 
M 2 D = rh 2 QL + m 2 D + m 2 Hl + A 2 b (43) 

and the initial condition for S, in the case of MSSM, is given by 

So = (m 2 H2 ~ m 2 Hi ) Q + Y,(m 2 Q + m 2 D + m 2 E -m 2 L - 2ml) Q ( 44 ) 

gen 

In the case of the universal boundary conditions the S contribution vanishes 
identically as required by the absence of gravitational mixed anomaly [p4fl . In this 
case S vanishes at any scale as can be deduced by integrating its corresponding 
renormalisation group equation, S(t) = ^^-Tr [Ym 2 ]. It is obvious from the 
above equations that if we start with the same or almost equal h t , h b couplings, 
the two Higgses are going to evolve in a similar way. As discussed in the case 
of Yukawas, again small differences of U(l) factors in the scalar mass system 
cannot be enough to generate a difference of m H , and m 2 H2 running masses at 
mw- Here, however, in contrast to the h t — h b case, one has two possibilities 
of generating m 2 Hi — m 2 H2 > m 2 z at rri\y- The first is due to the difference of 
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the Yukawa couplings ht t b- In fact the low tan/3 regime is the limiting case of 
this possibility. One uses for economy universal b.c.'s for scalar masses at Mjj, 
therefore m° Hi = m° H2 and the difference arising at raw is only due to the fact 
that h t ^> hb. Here, if hb ~ 0{h t ), universal conditions cannot possibly lead to 
m 2 Hl 7^ m 2 H2 , at least not for the biggest part of the parameter space. A complete 
analysis of the parameter space (mi/2, rn 3 / 2 , fi, ...) at Mu is obviously required 
to determine those regions which are compatible with the radiative symmetry 
breaking and the phenomenological expectations. The isolation of the RSB- 
compatible regions of (777.3/2, wii/2) will have a further consequence in low energy 
phenomenology: it will evidently constrain the region of the supersymmetric 
sparticle spectrum masses. Such a constrained situation may serve as a guide in 
the experimental searches for their detection. 



Up to know, analytic solutions exist only in the case where h t ^> /i{, iT ....|]16[ |17 . 
In order to solve the above equations in the case of h t ~ hb, we form the differential 
equations for two sums of the above system, namely for M\,}A 2 D . Here we will 
work only in the limiting case where A^t — ► 0. In fact, by recalling the same 
arguments used in the solution of scalar masses for the case h t 3> /i 6i0 []T7| we 



can conclude that in the IR-limit the contributions of the trilinear parameters 
At, Ab do not play important role in the final solutions for the scalar masses as long 
as mi/2 < fn 3 /2- Therefore, to simplify the subsequent analysis, we drop Au,Ad 
terms (the extension of the solution to the most general case is straightforward). 
It is easily then observed that one can write the equations for the sums of scalar 
masses in the following form: 

M^ + MM-G° m y (45) 
i^h 2 + 6M 2 D h 2 b - G° D m 2 } (46) 



dMl 


1 


dt 


8vr 2 


dM 2 D 


1 


dt 


8tt 2 



where Gfj = Gq + Gh 2 + G\j^ and G° D = Gq + Gh ± + Ggc In this final system, 
we also observe that the S contribution is no longer present in the equations. 
Indeed, the two sums formed above include the partners of trilinear terms of the 
superpotential. Due to the U(l) invariance of the Yukawa Lagrangian, the S 
term is identically zero in the sums. 

To simplify the above coupled equations we make the following transforma- 
tions 



Ml = tX,t = exp^Jlhjdt'} 
M 2 D = ay, a = exp{^ jlh 2 dt>) 



(47) 
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Then Eqs( 



can be written in the form 
dX 



T 



a 



dt 

dy 

dt 



1 da 

6~dl 
ldr 
Q~di 



y 

X 







8tt 2 


m l/2 






8tt 2 


™i/a 



(48) 
(49) 



Using the formulae of the previous section, the above equations can be written 
as follows 



d_ 
dt 



X 

y 



H(t) 



X 

y 



m 



S7T-' I 4 



1/2 



(50) 



with 



2/o ( 



17 
12 



Zo(£)» 







(51) 



We may use Green's functions techniques to solve this non-homogeneous system. 
For this we need first the solution of the homogeneous part. 
If we define the function h{u) as follows 



h(u) = 

the homogeneous part becomes 

dX = - 

dy- 



Xq 



1 



u 



17 
Hi 



(52) 



du 



lOv^T^ 
1 du 1 



10 VvF+uh(u 



h(u)y{u) 
X(u) 



(53) 
(54) 



There is a rather complicated dependence on u through h{u) , however, it is 
possible to use new parameters which may simplify the system. For example, if 
we define 

dQ = ~~^= (55) 



lOx/^T 



u 



and 



we can easily find that 



1 



VT+u 
1 



sin2(f) 



Q = -ln(tari(f)) 



(56) 
(57) 
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Now the homogeneous system becomes 

dX 



dQ 



h{Q)y{Q) (58) 



dQ h(Q) 



If we seek solutions of the form 

X = Exp{ J hAdQ), y = Exp{ J S-dQ) (60) 

then it can be shown that the unknown functions A, B satisfy the relation B = 
1/A while we find that A obeys the first order Differential Equation 

^ = M±-hA) (61) 
dQ hA 

The latter, is a first order Riccati type Equation of the general form A + A 2 — 
f(x) = whose solutions are known only for particular forms of the function f(x). 
In our case, x = j h(Q)dQ and f(x) = l/h(Q) 2 , which is a rather complicated 
function of the new parameter x. Therefore exact analytic solutions are possible 
only in particular limiting cases. We can study the above equation by geometric 
methods and obtain useful information about its solutions. Such a mathematical 
analysis, however, is beyond the purposes of this paper.Q We will present soon, 
two limiting cases which are of particular interest from the physics point of view. 
Before we proceed to the presentation of the interesting limiting cases let us 
exhibit some interesting properties of the equation (^). The above equation 
remains invariant when A is shifted by a function g(Q), i.e. A — >• A + g, which 
obeys the equation 



g(Q) 



g Exp(-2jAhdQ)dQ 

iQ 

(62) 



go J hExp(-2 J AhdQ)dQ 

go i 



1 + g J hX~ 2 dQ' X< 



Furthermore, the equation A(u) 2 — h(u) 2 = defines a curve which separates the 
areas where the derivative of A(u) changes sign. Therefore above (below) that 
curve, A(u) is an increasing (decreasing) function of the parameter u. This is 
depicted in figure 2. Let us also discuss the possibility of searching for a simplified 



solution of the system (|5g,[5l|). We observe from figure 2 that for approximately 



^ee for example 
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equal couplings, the function h(u) is approximately constant within the biggest 
range of the parameter u. In figure 3 we compare the numerical solution for the 
homogeneous part of the differential system with the one obtained by considering 
h(u) ~ const., assuming non-universal boundary conditions. We can conclude 
that the approximate analytic solution is in good agreement with the numerical 
one for ht/hj, < 1.1. In what follows, we will consider approximated formulae in 
the interesting case of approximately equal top — bottom Yukawas. In terms of 
the variable Q the complete system is written 






i 

HQ) 



HQ) 






(63) 



where where Gjj,d{Q) are properly defined through (|50|) . 
With the use of the diagonalizing matrix, 



U{h) 



V2~h 



h -h 
1 1 



(64) 



we can write the Differential Equations as follows 




1 





-1 
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1 1 
1 1 




(65) 



(where the new vector Q T = (X, y) is related to the old one Q$ = (X, y) with 
the matrix U, Q = U~ 1 Qq). For the case under investigation, for example, this 
form is particularly useful. In the limit h t fl ~ ^&,o, the parameter «> 1, thus the 
second term in the right hand side may be considered as a first order perturbation. 
There is another interesting limiting case which can be also easily treated, i.e., 
when i/u <C 17 in the range of integration. In that case (which corresponds to 
h t ^> ht), we may define the new variable P = ln(h) and the transformation 



V 



1 

V2 



(66) 



the differential system can be written as 
17 




1 




1 

1 




-17/-1 




(67) 



This particular case is treated more easily if one ignores coupling. Partic- 
ular expressions are found in the literature []T6|, 17|. Thus, in the present work 
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we restrict ourselves in the approximate bottom-top equality. It can be easily 
concluded in that case that u ^> 1 always and the function h(u) is approximately 
constant, h(u) m (|^) 6//5 - In this case Eq(|67D can be integrated, leading to the 
following result 

X \ „ „ „ I / Xo \ rQ ,„,„ r „ x , / G 



y 



ExplHoiQ-Qo)]^ [ y)+ J Qa dQ'Exp[-n (Q'-Q )] J ) | .(6 



where, due to the property H 2 = ( ^ ^ ] , the exponents can be expanded as 
follows 



tt \n-i in n M ut V " VO "O^W ~ VO ,„ n v 

Exp[H {Q - Qo)\ = ( (-1) | (69) 



ch(Q - Qo) h sh(Q - Q j 
h^ 1] sh(Q - Q ) ch{Q - Qo) 

In the large tanf3 regime and taking into account the Q(u) dependence of Eq(|5 
for the present purposes the above solutions can be approximated as follows 



M 2 r 11 

— T « O (^ + ^D)p+(^-^D)- +6/2</r> (70) 
m Z { P ) 

M 2 n a (V , _ _ , . ,1 



7T7- (^ + ^d)p-(^-^d)- +ei/2</ CT > (71) 

with p = ( fanto ^ 5 '' "*" n ^ e ^ as ^ ^ wo ex P ress i° ns we have used the convenient 
parametrisation .M \u,d) = i(u,D)m 2 and m^ /2 = £i/2m 2 Q , while = £h 2 +£q+& c 
and £d = £m + £q + Finally the integrals < J T)(7 > are given by 

< >= if w I* 1 + "< ( '»^ + - "(owo}* (72) 

and similarly for < I a > with the replacements r — > cr and — > Go, ^ — » ~ 
while z/ is given by 



2 



(73) 



VVi + M + 1/ 

The integrals < 7 T)CT > look complicated, but in fact they can easily be performed 



using definitions ([12] , 18) to express the functions p, u, r, a in terms of the integra- 



tion variable t'. As a matter of fact, for the case of interest (hf, ~ ht), u ^> 1 the 
integrals can be simplified further, as u — > 1. Then, instead of (72), we simply 
write 

< j t >^ r / ... ... at , < I„ >~ a / , ... at 74 

it T(t')p(t') ' A (J(t')p(t') K ' 
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Once the solutions of the sums are obtained, we may substitute back into the 
equations for the individual masses. A simple way to solve these equations is 
achieved if we formally integrate the sums to obtain 

Ml - M 2 U>0 - Cu{t)m\ /2 = -6Ju - J D (75) 
M 2 D -M 2 Dfl -C D {t)m\ /2 = -Ju-6J D (76) 

with J i{t) = J hj-Adj-dt, I = U,D. Now, the unknown integrals Ji(i) can be 
expressed in terms of the already calculated sums, their initial conditions and 
known gauge functions, from the simple algebraic system ( 75^) . Then, the 
higgs mass parameters for example can be given from 

m 2 Hl = (£ Hl + C H (t)^ /2 )m 2 Q - 3J D (t) + I s (77) 
m 2 H2 = (U + C H (t)^ 1/2 )m 2 - 3J v (t) - I s (78) 

with Is representing the integral of the ^-contribution in the case of non - uni- 
versality. Similar expressions hold also for the other scalar masses of the third 
generation. 

In order to examine the radiative breaking of the electroweak symmetry, we are 
finally interested for the mass parameters fi 2 = m 2 Hj + fi 2 . The initial value of fi is 
also unknown, but this arbitrariness can be eliminated by using the minimization 
conditions and determine its value in terms of known quantities. For the large f3 
we can write 

Tfl 2 I 

^ K -™«» + t^p ~ 2 m| (79) 
Notice that with the use of the relevant renormalization group equation for /i and 
(f79l), we may determine the initial condition for the /x parameter which may be 
useful in several cases. 

Let us conclude this section with some remarks concerning the advantage 
of the above semi-analytic procedure with respect to the standard numerical 
integration. First we have shown that the useful quantities of the minimal susy 
model (Yukawa solutions, scalar masses, tan/?, /i-term e.t.c.) can be expressed in 
terms only of known gauge functions and the parameter u. The value of the latter 
can be obtained easily at any scale in terms of the top and bottom couplings at 
the unification scale through the hypergeometric function defined in fln|) . In fact, 
in the limiting case we discuss in this work, the parameter u can be expressed 
by a simple, approximately linear relation in terms of its initial value, thus a 
lot of information about the low energy parameters can be easily obtained by a 
straightforward calculation of the involved parameters. Furthermore, the scalar 
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spectrum depends on two more new quantities, i.e. the integrals J T|(T , which can 
also be easily computed or fitted by simple mathematical functions for any given 

h(t,b)o set - 

A final comment we wish to make here concerns the possible effects of the 
At j, parameters in our solutions. We have already noted that, the A t — A^ renor- 
malisation group equations can be solved independently in terms of the gauge 
and Yukawa couplings. In fact, the A t b coupled differential system is exactly 
the same with that of the scalar sums ( |i6|) with the proper replacement of the 
coefficients Gu,d- After replacing their solutions A tj b = ft,b(At,o, ^4&,o) — /[ ff mi/ 2 , 
into the system of the scalar masses, one finds that the 8a, corrections are of the 



form (see also IE6 



5 Ai = S ix Al + 5 i2 A m 1/2 + 5 h m 2 1/2 (80) 
The coefficients S~ n „ involve the two more integrals similar to I r „ defined above. 

{1,2,-5) ) 

Remarkably, in the case of Yukawa couplings close to their infrared value, the 
coefficients multiplying the Aq terms are found to be small 5i 12 <C 1. This 
result has also been confirmed elsewhere ^6] , where all terms involving the initial 
condition A have been shown to be multiplied by a factor 1 — (-j^) 2 , with h{ xd 
the fixed point value of the top Yukawa coupling. For reasonable A , values, 
(i.e., not much larger than the gravitino mass), and staying close to the fixed 
point /if-value, the above corrections do not modify substantially the present 
analysis. 



22 



5 Results and Conclusions 



We are ready now to present our numerical examples using the obtained scalar 
mass formulae. Although the formulae of the previous section look rather compli- 
cated, in fact for most of the cases of interest, they can be simplified considerably. 
As a matter of fact, the only "new" quantities that enter the various scalar masses 
are Ju,d which mainly depend on the two integrals I T , l a and the initial bound- 
ary conditions For example, in the case of h t ~ hf, couplings close to their 
non-perturbative regime, the Ju,d quantities can be easily simplified as follows 



^~^{3.0 + e c + (CV(t)-/ T )6/ 2 } 
J D ™\{w + ez + {C D (t)-h)£ 1/2 } 



(82) 



where e^, eg represent small corrections of the order 0(~ 3%) which can be easily 
extracted from the previous formulae. Furthermore the Cu,d 
functions are easily calculated from the formulae 



EnCn' D gauge 



Cn(t) = £ 



- (a$(tx) - a?(t) 



n 



1,2,3 



Thus Cy ~ 7.1 + 6.7 + 0.5 = 14.3 and similarly for Cd- Finally, for the region 
of interest the two integrals I T a are easily calculated. For h to ^3.0 for example 
we find I T ~ 5.7, while as h^ Q < h to , we obtain I a > I T . In the case of equal 
bottom-tau couplings the two expressions (81], 82]) are almost identical (up to 
U(l) differences in Cu,d gauge quantities and the eg, eg parameters which might 
differ due to non-universal initial boundary conditions). In the limiting case of 
equal b.c's at the GUT scale these quantities are equal and our expressions reduce 
to those of similar investigations of references |26| , 

Let us now turn to the numerical investigation of the scalar masses. Using the 
the minimization conditions for the potential, it has been shown in section 2 that 
in order to achieve Radiative Symmetry Breaking at low energies, the constraint 
^ m Hi 2 = — m ^2 — m z should be satisfied. After some simple algebra we 
arrive at 



ml 



{2(frx-&&)+3(£«<'-&0} + 2is 



T 



+ 



5 
3 

To 



h h p 



Cd,c) + I(tt\ £l/2 



34) 
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with the I aT integral being the difference I a — I T . In the above formulae, all 
functions depending on the parameter u are finally expressed in terms of the 
initial value Uq, which depends only on the ratio of Yukawas, ht,o/hb,o at Mqut- 
Therefore, the constraint m 2 H — rn 2 H2 > m\ has been reduced down to a simple 
algebraic inequality of the form 

afffx + Kh 2 + lii/2 + 5>Zz (85) 

with £z — ( m z/ m o) 2 while a,/3, 7 and 5 can be easily read from flgl]) . (The 
other constraints can also be easily converted to simple equations of this latter 
form). Notice that in the limit hf, — > h t the integrals I TjCT are approximately equal 
while C u c,Cdc differ only in the U(l) factors. Therefore, in the case of non - 
universality, in general we get 7 < a, (3 from (|85|) , and the main dependence is on 
£hi,2 parameters. For h t > h b at M v we get I aT > 0, therefore the coefficient in 
front of £1/2 is always positive. Furthermore, for universal boundary conditions 
= £,h 2 ,£,u c = (,d c the first line in the RHS of equation (|84]) becomes zero. In 
this case a lower bound is put on value of £x/ 2 to satisfy (|85|) . In any case, this 
condition can be naturally satisfied for non - universal boundary conditions and 
in particular if £ Hl > C,h 2 ^u c > £,d c - Thus, once a particular point (£hi> £ff2> C1/2) 
is chosen, this condition may be seen as a lower bound on the scale 7713/2, i.e. 
m 3/2 > nT>z/(c^CHi + P£,h 2 + 7G/2 + S) 1 ^ 2 . In fig.^, this constraint is represented by 
a two dimensional surface as a function of the £#1 and £i/ 2 parameters choosing 
two specific values for £ H2 . Evidently, for small £i/ 2 the constraint is satisfied 
only for f Hl > 

We may further examine the Higgs mass term m 2 H2 which should be negative 
in order to break the SU(2) symmetry. If we substitute Jy from the previous 
equations in (|77|), (ignoring for simplicity small corrections of the order (9(2 — 
3%)), we arrive at the result: 

^ w ^{(17^ 2 + 3^ 1 )-3(5^ + 6^-^)} 

35 

- \{G Q + C u c - ^C H - JJ6/2 (86) 

After substituting the relevant functions by their numerical values, one can easily 
check under what conditions the required inequality (777^, < 0) is satisfied. For 
example assuming universal conditions for all scalars except 777^ , a lower bound 
is put in the parameter £i/ 2 only if £ Hl > 4 i.e., when m Hl > 2m 3 / 2 - 
Similarly, for the t c squark running mass we get 
fh 2 I 

-T ~ ^{(23^c + 2^)-2(5e Q + 6U-U)} 

/77-q oO 
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C u c — I T }£l/2 



(87) 



Now, if we denote with m e ^ v ' the lower experimental bound on the stop-mass, 
it is rather simple to investigate the constraint m\ c > {fh 2 c ) exp ' from the above 
formula. In fact, this can be also converted to a bound on the parameter 777.3/2 as 
in the case of Eq. fl85|) . In most of the parameter space however, this bound on 
m\ is lower than the corresponding obtained from 5m 2 Hi2 condition. In fact, the 
(mf c ) exp ' - bound excludes only regions with £y 2 ^ 1- 

Nevertheless, interesting bounds - in addition to ^ > - arise mainly 
from the m&c - mass, given by 



For many interesting regions of the £j parameter space the bounds obtained from 
this latter mass are non -trivial. Assume for example a moderate case where 
£q = = £# 2 = ~. (This would correspond to the relation 777,3 ~ 0.7m 3 / 2 for the 
relevant mass parameters). Then the above formula implies the following relation 
between £1/2, 6fi 



where m e b xp is the experimental lower bound of the sbottom mass. Now, if we 
imply a sensible value for (say m# a ~ 2777,3/2) so as to satisfy naturally the 



ratio mj£ p /rriQ is. This in turn implies the constraint > (0.55 — 0.65)m 3 /2. Of 
course this holds for the particular set of &'s chosen, but a detailed investigation 
can be easily done for the whole parameter space. 

We should finally point out that the above results do not include contributions 
from a large /i T -Yukawa coupling which in various models is also predicted to be 
of the same order with h t) b- A complete account of these effects in the described 
analytic procedure up to now was not possible. However, corrections of this 
Yukawa coupling are not expected to have a significant impact on the large part 
of the parameter space. Nevertheless, it will be a rather interesting question if 
the above analytic procedure can also be extended to include corrections of h T , 
in the case where ht ~ hb ~ h T . 

In conclusion, in this paper we have presented analytic expressions for the 
Higgs mass parameters m 2 Hl and squark masses of the third generation, when 




(88) 




(89) 



condition ((S3]), we should have at least £i/ 2 > (0.3 — 0.4), no matter how small the 
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both top and bottom Yukawa couplings are large, as predicted by many unified 
theories. In the derivation of these expressions we have used exact analytic solu- 
tions of the top-bottom (h t , hb) running Yukawa couplings ( expressed in terms 
of hypergeometric functions) which are involved in the evolution equations of the 
scalar masses. We have shown that, in the limit ht, hb 3> h T , ... the coupled differ- 
ential system of Higgses and scalars quarks is reduced down to a simple first order 
Riccati Differential equation. Specific numerical examples have been presented in 
the interesting case of large tan/3 ~ O(60), using values for the top coupling close 
to its non-perturbative regime. For these ranges of the parameters ht, tan/3, and a 
bottom mass Pd 4.25GeV, we predict a running top mass m t ~ (176 — 177)GeV 
corresponding to a physical mass mf h (184 — 185)GeV. These predictions can 
change substantially, if corrections from sparticle exchange loops (being propor- 
tional to tan/3) are taken into account. They depend crucially on the value of 
the \x parameter, while they are mainly controlled by a linear combination of the 
g\vn- g and h 2 A t quantities, where nig, A t are the gaugino mass and the trilinear 
coupling respectively. It is possible of course to assume proper boundary condi- 
tions to suppress those corrections down to a (5 — 6) % of the bottom mass, but 
in principle they can be as large as 40% is specific regions of the parameter space. 
A final answer to this question can be given only when the underlying theory is 
known to pick up a particular point of the parameter space. For the time being 
we can reverse the procedure and constrain the arbitrary parameters using the 
experimentally determined range of the top mass. 

We have further presented simple expressions for the Higgses and scalars and 
discussed the constraints put by the radiative symmetry breaking mechanism, as 
well as, by charge and color protection on the initial conditions of the scalars. 
Thus the condition that the one Higgs (mass) 2 turns negative at low energies can 
be easily satisfied for natural values of the scalar masses as can be concluded from 
(|56j). Furthermore, the second Higgs is protected from large negative corrections 
at low mass scales, either if we impose m 2 H > m 2 H2 at the unification point, or if 
TO1/2 > m 3 / 2 . It is evident from ( [34]) that the specific boundary conditions on the 
squark masses mt^,m^, will play a role in the exact determination of the ranges 
of the above mass parameters. Finally mQ,m t c,rhbc squark masses also impose 
additional constraints on the (777,3/2, my 2) mass parameters. 
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FIGURE CAPTIONS 



Figure 1. Plots of the ratio 773 = (tan/3 /tanfio) 2 ^ vs the parameter u for 
selected ratios of the Yukawa couplings htfi/h^Q. 

Figure 2. The curve A(u) = h(u) separating the two regions with A > and 

A<o. 

Figure 3. Comparison of analytic approximate (upper curves) and numerical 
solutions (lower curves) of the ratios X(u)/ml, y(u)/m 2 0) for two cases of Yukawa 
couplings. (For X,y definitions see section 4). 

Figure 4- Surfaces representing the lower tUq bound, in the parameter space 
£1/2 = m i/2/ m oi £#1 = m Hi/ m o f° r £h 2 — 0.64 (lower) and — 1 (upper case). 
If £,h 2 > £hd surface points exist only for £y 2 ^> 1. 
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Figure 4 (from different viewpoints) 
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